Neutrino emission from spin waves in neutron 
spin-triplet superfluid 



L. B. Leinson 

Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation RAS, 
142190 Troitsk, Moscow Region, Russia 



Abstract 

The linear response of a neutron spin-triplet superfluid onto external weak axial- 
vector field is studied for the case of 3 i-2 pairing with a projection of the total angular 
momentum rrij = 0. The problem is considered in BCS approximation discarding 
Fermi-liquid effects. The anomalous axial- vector vertices of neutron quasiparticles 
possess singularities at some frequencies which specify existence of undamped spin- 
density waves in the Cooper condensate. The spin waves are of a low excitation 
energy and are kinematically able to decay into neutrino pairs through neutral weak 
currents. The calculation predicts significant energy losses from within a neutron 
star at lowest temperatures when all other mechanisms of neutrino emission are 
killed by the neutron and proton superfluidity. 
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When the temperature inside a neutron star core has dropped below the crit- 
ical temperature T c for spin-triplet neutron pairing, the bulk baryon matter 
is expected to develop a superfluid condensate of neutrons p]- [6], which has 
thermal excitations in the form of broken Cooper pairs. It is generally accepted 
that neutrino emission from the pair-recombination processes in the neutron 
triplet superfluid dominates the neutrino emissivities in many cases [?]. Ac- 
cording to the minimal cooling paradigm [8], [9], along with lowering of the 
temperature the star continues to lose its energy by radiating low-energy neu- 
trinos via the recombination processes untill enters a photon-cooling epoch at 
T ~ 0.1T C . At this latest stage of the neutrino-cooling, all other mechanisms 
of the neutrino emission from the inner core are suppressed greatly by the 
neutron and proton superfluidity [10] . 



In this paper we calculate the new neutrino emission process that can dom- 
inate in the temperature range typical for the final stage of neutrino cool- 
ing. Namely, we consider the weak decays of spin density excitations in the 
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superfluid triplet condensate of neutrons. Previously spin modes have been 
thoroughly studied in the p-wave superfluid liquid 3 He [16]- [20]. The pairing 
interaction in 3 He is invariant with respect to rotation of spin and orbital 
coordinates separately. This admits spin fluctuations independent of the or- 
bital coordinates. In contrast, the spin-triplet neutron condensate arises in 
the high-density neutron matter due mostly to spin-orbit interactions which 
do not possess the above symmetry. Therefore the results obtained for liquid 
3 He can not be applied directly to the spin-triplet neutron superfluid. 

Recently spin waves with the excitation energy smaller than the superfluid 
energy gap was predicted to exist in the superfluid spin-triplet condensate 
of neutrons [21 J. The neutrino decay of such spin waves could be important 
for thermally-emitting neutron stars, presumably cooling through the com- 
bination of neutrino emission from the interior and photon cooling from the 
surface, the latter is responsible for their observed thermal emissions [22]. 

Let us remind briefly the theory of spin density excitations in the condensate. 
The triplet order parameter, D = D a/3 (n), in the neutron superfluid represents 
a symmetric matrix in spin space (a, (3 =t, \) which can be written as D (n) = 
Ab&g, where <x = &2, 03) are Pauli spin matrices, and g = i&2- The gap 
amplitude A is a complex constant (on the Fermi surface), and b (n) is a real 
vector which we normalize by the condition^ 

/ ^ 2 (n) = 1 . (1) 
J Air 

The angular dependence of the order parameter is represented by the unit 
vector n = p/p which defines the polar angles (8, tp) on the Fermi surface. 

The neutron pairing at high densities involves a mixing of 3 P2 and 3 F 2 channels 
Conventionally this mixing is not taken into account in estimates of 
neutrino emission from the neutron superfluid. Accordingly, throughout this 
paper we consider the case of 3 P2 neutron pairing, when quasiparticles pair 
in the most attractive channel with spin, orbital and total angular momenta, 
s — 1,1 — l,j — 2, respectively. Then the pairing interaction, in the most 
attractive channel, can be written as pQ 

"Tap,* (P, P') = -V ip,p') £ (b m »*$) ($* b ^( n ')L ' ( 2 ) 

m,j 

where V(p,p') is the corresponding interaction amplitude; v = ppM*/7T 2 is 
the density of states near the Fermi surface, and h mj (n) are the vectors in 



In what follows we use the Standard Model of weak interactions, the system of 
units h = c = 1 and the Boltzmann constant ks = 1- 
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spin space which generate standard spin-angle matrices so, that 

b mj (n)ag = Y ( ^ a ^\ lm s) (Hrn s mi\2mj) Y Xjmi (n) . (3) 



m a +mi=mj 

These are given by 



b = Y / T/2(-ni,-n 2 ,2n 3 ) , bi = -y / 3/4 (ra 3 , m 3 , rii + m 2 ) , 

b 2 = ^3/1 (m + m 2 , im - n 2 , 0) , b_ m . = (-) m > b^. , (4) 

where n\ = sin # cos ip, n 2 = sin # sin ip, n 3 = cos#. The vectors are mutually 
orthogonal and are normalized by the condition 



(in 
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bm'.brrtj — S mjm '. (5) 



We will focus on the condensation with rrij = which is conventionally con- 
sidered as the preferable one in the bulk matter of neutron stars [I], [8], [9]. 
In this case b (n) = bo (n). 

Making use of the adopted graphical notation for the ordinary and anomalous 
propagators, G = — , G~(p) = — « — , = » « , and F^> = - - , it 
is convenient to employ the Matsubara calculation technique for the system in 
thermal equilibrium. Then the analytic form of the propagators is as follows 



G(p m ,p) = G(p m ,p)5 al3 , G (p m ,p) = G (p m ,p)6 a/3 , 
F {1) {Pm, P) = F (Pm, P) bag , {p m , p) = F {p m , p) gab , (6) 

where the scalar Green's functions are of the form G~ (p m , p) = G (— Pm , — p) 
and 

GiPrniP) = J" ^ , F( Pm ,p) = . (7) 

fm ' p "m ~ p 

Here p m = in (2m + 1) T with m = 0, ±1,±2... is the Matsubara's fermion 
frequency, and e p = p 2 / (2M*) - p%/ (2M*) with M* = p F /V F being the ef- 
fective mass of a quasiparticle. The quasiparticle energy is given by E 2 = e 2 v + 
A 2 b 2 (n), where the (temperature-dependent) energy gap, Ab (n), is anisotropic. 
In the absence of external fields, the gap amplitude A is real. 

Finally we introduce the following notation used below. We designate as 
(to, n, q; T) the analytical continuations onto the upper-half plane of com- 
plex variable u of the following Matsubara sums: 

i X x> (u n , n ) q; T ) = T Yl 2 / d£pX ( Pm + Wn ' P+ f) X ' ( Pm ' P ~l) ' ^ 
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where X, X' e G,F,G , and u n = 2mTn with n = 0, ±1, ±2. ...These are 
functions of u, q and the direction of the quasiparticle momentum p = pn. 

We will focus on the processes with u> 2 < 2A 2 b 2 and with a time-like mo- 
mentum transfer, q 2 < u 2 . In this case the key role in the response the- 
ory belongs to the loop integral Xpp. A straightforward calculation yields 
1 FF (n, u, qn; T) = X (n,w; T) + O (q 2 V£/uj 2 ), where 

/•oo ^£ A 2 F] 

X (n,w;T)=2/ — — — — — 1 tanh — , (9) 

Jo E AE 2 - (oj + i0) 2T 

with 

E = ^Je 2 + A 2 b 2 . 

Insofar as q 2 Vp/u 2 1 and q 2 Vp/A 2 C 1 we will neglect everywhere small 
corrections caused by a finite value of space momentum q. 

We are interested in the linear medium response to the external axial-vector 
field. The field interaction with a superfluid should be described with the aid 
of four effective three-point vertices. There are two ordinary effective vertices 
corresponding to creation of a particle and a hole by the field (These differ 
by direction of fermion lines), and two anomalous vertices corresponding to 
creation of two particles or two holes. 

The anomalous effective vertices are given by the infinite sums of the dia- 
grams taking account of the pairing interaction in the ladder approximation. 
In general, the ordinary vertices are to incorporate also particle-hole inter- 
actions. Since Landau parameters for the particle-hole interactions in asym- 
metric nuclear matter are unknown, we simply neglect the Fermi-liquid effects 
and consider the effective vertices and the pair correlation function in the BCS 
approximation. In this case the ordinary axial-vector vertices of a particle and 
a hole are to be taken as <r and <r T , respectively. 

Given by the sum of ladder-type diagrams [25] , the anomalous vertices are to 
satisfy the Dyson's equations symbolically depicted by graphs in Fig. CD 

The vertex equations are to be complemented by the gap equation which is of 
the following form [lj in the case rrij =0: 

A (p) = -— / dp'p' 2 V {p,p') A Op') I ^-^1 tanh . (10) 



We are interested in the processes occuring in a vicinity of the Fermi surface. 
To get rid of the integration over the regions far from the Fermi surface we 
renormalize the interaction as suggested in Ref. [26]: we define 

V« (p,p'; T) = V (p,pO - V (p,p') (GGr) (p,p'; T) , (11) 
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>f^< y^^y 



Fig. 1. Dyson's equations for the anomalous vertices. The shaded rectangle repre- 
sents the pairing interaction. 

where the loop {GG~) n is evaluated in the normal (nonsuperfluid) state. In 
terms of the gap equation can be reduced to the following simple form 



l = _yM f^b 2 (n)A(n 

J 4:71 



(12) 



assuming that in the narrow vicinity of the Fermi surface the smooth functions 
V (p,p') and A (jf) may be replaced with constants. 

The function A (n) arises due to the renormalization procedure. It is given by 

1 



A (n) = - [°° de [ - 

2 Jo We 2 + A 2 b 2 



, Ve 2 + A 2 6 2 1 , £ N 

tanh — tanh — — 

2T £ 2T 



(13) 



Details of the analytic calculation can be found in Ref. [21] . where after the 
proper renormalization the equations for the axial-vector anomalous vertices 
are obtained in the following analytic form: 



Tf ) (n) = V (r) £ *b mj (n)g I ^- [Z GG -Tr (g (ah*) f, 1 



da! 



8tt 



=.(1) 



X FF Tr ((*b^) (*b) ^T/ 2) («xb)) 



— Z FF 2i (b^xb 



(14) 



T/ 2) (n 



= K^^(n) / ^ [z G - G Tr ((<xb m ,) gf t 

- XppTi ((a-b m .) («rb) rf ^ (as)) 

-—l FF 2i [b mj xb 



(2) 



(15) 
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The consistent solution to equations (fl2|) . (JHJ), and (fT5|) is found to be [21] : 
f W (n;w, T) — f (u, T) (e (frbj) + e* (*b_0) £ , (16) 

T< a > (n;w, T) = f (u, T) g (e (<xbO + e* (d-6_0) , (17) 

where e = (1, — i, 0) is a constant complex vector in spin space; and the func- 
tion f (oj,T) is given by the expression 

/ (w, T) = —1— t/f 4- / — b 2 X FF (b 2 , u>, T) (18) 

with 



X (u,T) = J^ (h\h x - b 2 ) A (b 2 ^T 



, dn 



4A 1 



(b*b0- (Kb)(bbi) i FF (b 



(19) 



In this expression the vectors b 1; b^, and b = b are given by Eqs. (jlj). As is 
well known, poles of the vertex function correspond to collective eigen-modes of 
the system. Eigen frequencies, u = u s , of such oscillations satisfy the equation 
X (w s ,T) = 0. Since we consider the time-like domain, q < u s , the correction 
caused by finite value of the wave number q is proportional to q 2 V F <C u 2 (0). 
We will neglect this small (positive) correction to the dispersion of spin waves. 



Thus the dispersion equation for the collective oscillations is of the form: 



A(b 2 ,u } T 



+2 



dn 
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bjbi)- (bib) (b^) I FF (P,u,T) =0 . (20) 



Before proceeding to the detailed solution of Eq. (1201 . let us note that the 
equilibrium order parameter is specified in this equation solely by means of 
the real vector b. Therefore this equation allows to obtain the eigen energy 
of a similar excitation in a superfluid 3 He-B (Balian-Werthamer phase). In 
this case the order parameter matrix is given by D (n) = Ab<xg with b = n, 
and the energy gap is isotropic [27]. The latter means that, in Eq. (1201) . the 
functions A and X FF are isotropic and can be moved beyond the integrals. 
Using also Eqs. (JTJ), (JH]) we find 
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thus obtaining the dispersion equation 



r dn 
J 4^ 



4A 2 



(bjbx)- (bln)(b,n) 



(22) 



which has the solution uj = y8/5A, in accordance with the energy of the spin 
wave in 3 He-B obtained in |20|. 



We now return to examination of spin waves in the spin-triplet neutron su- 
perfluid. A simple estimate of the excitation energy of the wave can be made 
using the angle-averaged energy gap in the quasiparticle energy, (^A 2 b 2 ^ = A 2 . 
In this approximation the functions A and Xpp are isotropic and can be moved 
beyond the integrals. Using also Eqs. (l2Tj) we obtain the equation 



dn 

Air 



4A 2 



(bJbO- (Kb) (^b 



(23) 



with b = b . This equation has the solution uj = A/\/5. This simple estimate 
shows that the energy of the spin wave excitation is smaller than the energy 
gap in the quasiparticle spectrum. Neutrino decays of such spin waves could 
be important for thermally-emitting neutron stars. Therefore let us examine 
the wave excitation energy more accurately. 

Since the functions A (b 2 ,uj } T^j and Iff (& 2 >w>^) are axially symmetric one 
can integrate Eq. (|20|) over azimuth angle. We then obtain the equation 



J dn 3 (l - 3n 2 ) A (n 3 ) 

= -3^3(^(1 + ^)- n 2 (l-n 2 ))z K^,T) • (24) 

It is convenient to define the dimensionless variables y = A (T) jT and Q = 
uj I (2A). We can get a clear idea of the behavior of the left side of Eq. ( 124|) 
as a function of Q near the critical temperature. Obviously, y tends to zero 
as the temperature nears the critical value, T — > T c . In this case we find 
A (%; T) — > —>ty 2 b 2 , where 

x = — j°° % ftanha; - x (l - tanh 2 x)) = 5.32 x 10 -2 , 
16 jo x ^ v ' ' 

In the same limit T — > T c , the function X , tends to 

ov d ' ' y; 8 V6 2 - n 2 
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We obtain the dispersion equation of the form 
2(l + n§) 



fit 



dn? 



1 + 3n 3 



2W S 



o 



i n| (1 - n§) 



1 + 3n 3 



2Qi 



-0.02?/ 



According to the above estimate fijCl. Therefore the equation can be solved 
by iterations. To the lowest accuracy we find 

Q 2 S ~ 4.28 x 10" 2 - 1.76 x lQ- 2 y 



For arbitrary y > the dispersion equation requires numerical computations. 
In Fig. [2j the dimensionless frequency Q s = u s (0) / (2A) of the collective spin 
oscillations is shown versus the dimensionless parameter y = A (T) /T. 




Fig. 2. Dimensionless frequency £l s = oj s (0) /2A and function A vs. the temperature 
parameter y = A (T) /T. 



Just below the critical temperature, y — >■ 0, we obtain Q s ~ 0.21, or, identi- 
cally, u s (T — > T c ) ~ 0.42A. The spin wave frequency slightly decreases along 
with lowering of the temperature and tends to a plateau at Q s ~ 0.15, what 
is equivalent to u s (T — > 0) ~ 0.3A. 

Since the spin excitations occur in the time-like domain of the transferred 
energy and momentum the decay into neutrino-pairs is kinematically allowed. 
We examine the neutrino energy losses in the standard model of weak in- 
teractions. Then after integration over the phase volume of freely escaping 
neutrinos and antineutrinos the total energy which is emitted per unit volume 
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and time can be obtained in the form (see details, e.g., in Ref. [28]): 

6 " 192^ X i g exp (u/T) - 1 im Uwcak K qj V*"*" * ^ 

(25) 

where Gp is the Fermi coupling constant, M v = 3 is the number of neutrino 
flavors, = (w, q), (a?) is Heaviside step function, and n^ ak is the retarded 
weak polarization tensor of the medium. The latter can be found [21] using 
the linear correction to the Green function of a quasiparticle G' = SG/SA 
caused by the external field A. In the BCS approximation, it is given by the 
diagrams shown in Fig. [31 and can be written analytically as 



G* = 



Fig. 3. Correction to the ordinary propagator of a quasiparticle in external field. 

G' = GG & + FF (<rb) <r (<xb) + GF f {1) g (<xb) + FG (<xb) gt (2) , (26) 
where GG = G (p m + u n , p + q/2) G (p m , p - q/2), ect. 

In the axial channel, the pair correlation function can be found as the analytic 
continuation of the following Matsubara sum 

U i l(u n ,q)=T^J^Tr(a l G' J ) , (27) 

and the imaginary part of the weak polarization tensor can be written as 
ImIT(^ ak = 5 tJ,l 5 u iC\ Imll^, where Ca — 1-26 is the axial- vector weak coupling 
constant for neutrons. 

To the lowest accuracy in Vp 1, we may evaluate the polarization tensor in 
the limit q = 0. Then using Eqs. f[T6"j) . f lT?]) . f l2"B"|) . and f[2"T|) we find after some 
algebraic manipulations (Dependence on n and u in the integrand is omitted 
for brevity): 



TVi (oo)=A P J^ X - (X GG - b 2 X FF ) b l3 + b 2 l FF 



bibj 
b 2 



+ (<% - 5 i3 5 j3 ) f H b 1 FF 



(28) 



Jq=0 



Using the identity [Igg — tfTpF) = — 26 2 X , which can be verified by a 

V / q— >0 
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straightforward calculation, we find 

dn ( bjhj 3 . > \/6 



Imaginary part of this function consists of two contributions. The first one 
arises from the imaginary part of Xq (uS) and exists at u > 2b A. This part is 
responsible for the emission of neutrino pairs from breaking and formation of 
the Cooper pairs at thermal equilibrium and is already discussed in Ref. |21j . 
We are interested in the neutrino emission from the decay of spin waves. The 
corresponding contribution into the imaginary part of the axial polarization 
tensor (1251) arises from the pole part of the function / (co) at x ( w ) — 0, and can 
be found with the aid of Sokhotsky's formula, (% + iO) 1 = V (1/x) ~~ 27r£ (x)- 
Taking into account Eq. (120]) we obtain: 

Im / ( u ) = - W5£ [J] dn, (1 + 3,4) 1 ) J^J^-S - .1) (30) 

We remind that Xq is real function in the domain < u> < 2Ab, where the 
undamped spin waves exist. To simplify Eq. ( 130]) we use the fact that u 2 <C 
4E 2 . This allows to neglect u 2 in Eq. (J9j, and substitute 

t 1 f°° du * i yi> I 2 i 1 

In — / tanh — vr + 1 

26 2 io ( M 2 + 1 )3/2 2 

in Eqs. ( l29l) and ( 1301) . In obtaining this expression we used the change £ = uAb. 
Then we find 

Im / (w) = -tt-= " J (w - w s ) , 

V6 Jo dn 3 (l + n§)X 

and the imaginary part of the weak polarization tensor can be written as 

37T 

Im lC ak = - W (5 V - 5 t3 5, 3 ) -C\pu 5 (co - u s ) A (y) , (31) 

where 

(fidn 3 (l + 3nj)X " 



A(y) 



As shown in Fig. 2, this function increases smoothly along with lowering of 
the temperature and tends to plateau at A = 1.35 

Inserting Eq. ( l3~Tj) into Eq. (125]) and performing trivial integrations we find 
the neutrino emissivity due to spin wave decays (SWD): 

e SWD = ^G 2 f C\KpfM*T 7 F (y) , (32) 
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Here y = A (T) /T, and 

y 7 A (y) nj (y) 



exp (2y^ s (j/)) - 1 ' 
For a practical usage we reduce Eq. (132]) to the traditional form 

eswD . 6 . 3 x (£) (£) rJCiF (A <T) /T) , (33) 

where M and M* are the effective and bare nucleon masses, respectively; c is 
speed of light, and T 9 = Tj (10 9 K). 

The function F (y) depends actually on the only parameter y and can be easily 
evaluated making use of the analytic fits which relate fl s and A to y at any 
V > 0: 

0.2172 - 0.005% + 0.011 V + 0.0026?/ 3 



V s (y) 



a (v) 



1 + 0.053% + 0.0710j/ 2 + 0.0175?/ 3 
1.10382/ + 0.2589j/ 2 + 0.1063?/ 3 



1 + 0.5653?/ + 0.209V + 0.0780?/ 3 ' 
The maximum fit error is about 0.1% both for Q s (y) and for A (y). 

It is convenient to fit also the function y = A (T) /T for practical computa- 
tions. For this purpose we can adjust the simple expression suggested in Ref. 
[7], where 

„ ( T)s VT^ (0.7893 (34) 

and r = T/T c . In Ref. [7], the gap amplitude A (T) is defined by the relation 
A 2 = A 2 (1 + 3 cos 2 6), while our definition is A 2 = \A 2 (1 + 3 cos 2 9), i.e. the 
gap amplitude A (T) is \f2 times larger than the gap amplitude A (T) used 
in Ref. [TJ, and thus y (r) = (r). 

The result of numerical computation is shown in Fig. HJ where we compare 
the neutrino emissivity due to spin wave decays with the self-consistent emis- 
sivity from the pair breaking and formation processes (PBF), as calculated in 
Ref. [21]. For a comparison we demonstrate the PBF energy losses ignoring 
the anomalous weak interactions, as suggested in Ref. [7]. The temperature 
dependence of the emissivity is evaluated at p = 2 x 10 14 g cm~ 3 . We set the 
effective nucleon masses M* = 0.7 M, for simplicity. The critical temperature 
for neutron pairing is chosen to be T c = 10 9 i^. Total of the energy losses, 
e = cpbf + cswd, is also shown. One can see that the decay of spin waves 
into neutrino pairs is very effective at low temperatures, when other known 
mechanisms of neutrino energy losses in the bulk neutron matter are strongly 
suppressed by superfluidity. 

Let us summarize our results. We have studied the response functions of the 
triplet neutron superfluid onto external axial-vector field with taking into ac- 
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YKL 

PBF 

SWD / \ 




7.75 8 8.25 8.5 8.75 



igT [K] 

Fig. 4. Temperature dependence of the neutrino emissivity due to Cooper pairing 
of neutrons and due to decay of spin waves at p = 2 x 10 14 g cm~ 3 and T c = 10 9 K. 
For a comparison, we show the emissivity as obtained in Ref. [7] (dash line). 



count of anomalous interactions. The calculations are made in the BCS ap- 
proximation for the case of 3 P 2 superfluid condensate with rrij = 0. 

Our theoretical analysis predicts the existence of spin-density waves of a small 
excitation energy u s ~ (0.3 -j- 0.4) A in superfluid condensate. We have cal- 
culated the neutrino energy losses caused by the decay of the spin waves 
through neutral weak currents. Neutrino energy losses due to spin wave de- 
cays are given by Eq. (1321) . Because of a rather small excitation energy the 
decay of spin waves leads to a significant neutrino emission at lowest tem- 
peratures T <C T c , when all other mechanisms of the neutrino energy losses 
are killed by a superfluidity. Since the neutron 3 P 2 pairing occurs in the core 
which contains more than 90 percents of the neutron star volume, the decay 
of spin waves can affect the minimal cooling scenario. In the enhanced cooling 
scenario, the effect of spin wave decays could be also quite noticeable if the 
proton (spin-singlet) superfluid suppresses greatly the direct Urea process. 
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